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ABSTRACT 

We present a study exploring the impact of a starburst on the properties of the surrounding circum- 
galactic medium(CGM): gas located beyond the galaxy's stellar body and extending out to the virial 
radius(~200 kpc). We obtained ultraviolet spectroscopic data from the Cosmic Origin Spectro- 
graph(COS) probing the CGM of 20 low-redshift foreground galaxies using background QSOs. Our 
sample consists of starburst and control galaxies. The latter comprises normal star- forming and pas- 
sive galaxies with similar stellar masses and impact parameters as the starbursts. We used optical 
spectra from the Sloan Digital Sky Survey(SDSS) to estimate the properties of the starbursts, inferring 
average ages of ~200 Myrs and burst fractions involving ~10% of their stellar mass. The COS data 
reveal highly ionized gas traced by C IV in 80%(4/5) of the starburst and in 17%(2/12) of the control 
sample. The two control galaxies with C IV absorbers differed from the four starbursts in showing 
multiple low-ionization transitions and strong saturated Lyman-alpha lines. They therefore appear to 
be physically different systems. We show that the C IV absorbers in the starburst CGM represent a 
significant baryon repository. The high detection rate of this highly ionized material in the starbursts 
suggests that starburst-driven winds can affect the CGM out to radii as large as 200kpc. This is 
plausible given the inferred properties of the starbursts and the known properties of starburst-driven 
winds. This would represent the first direct observational evidence of local starbursts impacting the 
bulk of their gaseous halos, and as such provides new evidence of the importance of this kind of 
feedback in the evolution of galaxies. 

Subject headings: galaxies: halos — galaxies: starbursts — galaxies: ISM — quasars: absorption lines 



1. INTRODUCTION 

The evolution of galaxies is regulated by how and when 
they accrete gas, produce stars and supermassive black 
holes, and by how different feedback processes may regu- 
late this cycle. Similarly, the evolution of the intergalac- 
tic medium (IGM) will be regulated by these same feed- 
back processes which can photoionize, collisionally heat, 
and chemically enrich the IGM. The interface between 
the region where star formation has taken place and 
the IGM is the galaxy halo, or circum-galactic medium 
(CGM). We take this to be the region extending from 
the main stellar body of the galaxy outward to the virial 
radius at a few hundred kpc. It is through the CGM 
that the two-way communication between galaxies and 
the IGM occurs. 

One of the major ways in which feedback occurs is 
via outflows of gas driven from strongly star-forming 
galaxies by the en ergy and /or moment um injected by 
massive stars (e.g. IVeilleux et al.l I2QQ5L and references 
therein). These outflows are u b iquitous in high r e dshift 
galaxies (fShaplev et all 12003b lAdelberger et al.l 12005b 
IWeiner e t al. 20Q9h. At low redshifts the conditions 
needed to drive these large-scale galactic winds seem 
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to be present only in starburst ga l axies ([Heckman et al.l 
19901: iLehnert fc Heckmanl 119961 : IVeilleux et al.l 12005b 
Martini [20051 : iRupke et aLll2005l and references therein). 
However, while less common, these low redshift outflows 
can be studied in far greater detail than their high-z 
counterparts. 

While it is clear that intensely star-forming galaxies at 
all redshifts drive outflows, the effect of these outflows on 
the evolution of galaxies and the IGM has not yet been 
well-quantified. Are they sufficiently strong to propagate 
out into the CGM and IGM? In so-doing can they eject 
baryons from low-mass dark matter halos (dwarf galax- 
ies)? Can they shut down the delivery of gas through 
the CGM in more massive systems, leading to the ces- 
sation of star-formation and the migration of galaxies 
from the blue star-forming sequence onto the red quies- 
cent sequence? What role do they play in establishing 
the galactic characteristic stellar mass of ~ 10 10 ' 5 M 
which separates the low-mass mainly-star-forming and 
high-mass mainly-quiescent galaxy populations? 

Unfortunately, the investigations of outflows in local 
starbursts cannot directly address these questions be- 
cause they have been limited to regions in the inner most 
part of the CGM (corresponding to radii of a few tens- 
of-kpc or less). Moreover, the outflow speeds that are 
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directly measured in the molecular, atomic, and warm 
ionized gas are typically comparable to or less than the 
galaxy escape velocity, so that it is not clear whether this 
mate r ial is truly a wind or a fountain (e.g. IChen et al.l 
[201(1 iGrimes et al.ll2009l ). The very high specific ener- 
gies implied by the high temperatures of the X-ray emit- 
ting gas suggest that this gas is likely to e scape at least 
from the less massive starbursting galaxies ([Grimes et al.l 
l2005h .but again these measurements pertain to locations 
far inside the virial radius. 

Direct imaging of the bulk of the CGM to probe the ef- 
fects of winds is difficult at high-z (but see ISteidel et al.l 
2011), and is impossible with existing facilities at low- 
redshift. The best available probe of the properties 
of the CGM comes from QSO absorption-line systems 
(QALS), which can detect this tenuous material with 
great sensitivity (albeit along single sight-lines). Various 
authors have suggested the co nnection between Q ALS 
and starburst driven winds. IZibetti et al.l ([20071 ) ar- 
gued that models of metal-enriched outflows from star- 
forming/bursting galax ies could exp l ain th e origin of 
strong Mg II absorbers. iBouche et al.l (|2007[ ) found that 
2/3 of their strong Mg II absorbers (with rest-frame 
equivalent width Wa2796 > 2A) at z ~ 1 had starburst 
galaxies (with SFR=l-20 M yr -1 ) around them . This 
result was also later confirmed by iNestor et al.l (|201lh 
with the discovery of bright, starburst /post starburst 
galaxies in the field around 2 "ultrastrong" Mg II QALs. 
Another indication of this c onnection comes from the 
study bv lMenard et al.l (|2Qllh . where they found a strong 
correlation between the equivalent width of the Mg II ab- 
sorbers and the [Oil] emission-line luminosi ty of the host 
galaxy. In a slightly different experiment, ISteidel et al.l 
(2010) have used more distant galaxies located behind a 
sample of intensely star-forming galaxies at z ~ 2 to 3 to 
trace outflowing gas in absorption to radii of ~100 kpc. 

With the i nstallation of the C osmic Origins Spectro- 
graph (COS; iGreen et al.l ([20121 )) on the Hubble Space 
Telescope (HST) it has now become feasible to probe the 
CGM of low-redshift galaxies via QALS. Most of these 
transitions lie in the ultraviolet and the superb sensitivity 
of COS allows us to directly detect and characterize the 
CGM in a representative samples of lo w-redshift galax- 
ies tha t can be studied in detail (e.g. iTumlinson et al.l 
I2011b[ ). In this paper we utilize COS to directly probe 
the CGM of starburst galaxies and to learn the effect 
of starburst-driven winds on the circumgalactic gas. To 
that end, we have selected a small sample of galaxies from 
the SDSS that have undergone a significant starburst 
event within the last few xlO 8 year, and we have also 
selected an appropriate control sample of typical star- 
forming and quiescent galaxies. For all these galaxies we 
have obtained ultraviolet (UV) spectroscopic data with 
the COS, which can directly detect gas in the CGM via 
absorption lines seen in the spectra of background UV- 
bright QSOs 

Detailed descriptions of our sample selection criterion, 
the COS observations, and data reduction are presented 
in Section 2. The results and their implications are dis- 
cussed in Section 3 and 4. And finally, we summarize our 
findings and conclude in Section 5. The cosmological pa- 
rameters used in this study are Hq = 70 km s _1 Mpc -1 , 
VL m = 0.3, and Q\ = 0.7. Atomic data used throughout 



our anal ysis w as obtained from the published work by 
lMortod|2003h . 

2. OBSERVATIONS 
2.1. Sample 

Our sample was selected in a two-phase process. The 
first step was to identify those galaxies in t he SDSS Data 
Relea se 7 (DR7) "Main galaxy sample ([Strauss et al.l 
120021 ) that are under-going or have recently undergone 
a strong sta rburst. To do so w e followed the procedure 
described in IWild et "all ([20071 ) . They used a Principal 
Component Analysis (PCA) approach to improve spec- 
troscopic diagnostics for bursts of star formation for a 
bulge dominated sample. In this analysis, SDSS spectra 
from rest-frame 3750-4150A were decomposed in multi- 
ple eigenspectra and the resulting principal component 
amplitudes represent the amount of each eigenspectrum 
present in the galaxy spectrum. The first principal com- 
ponent (PCI) provides t he same informatio n as the tra- 
ditional D n (4000) index ([Balogh et al.lll999l ). It is a good 
proxy for the specific star formation rate of the galaxy, 
averaged over the past few Gyr. The second principal 
component (PC2) provides additional information on any 
excess (or deficit) of B aimer absorption, over that ex- 
pected for the 4000 A break strength of the galaxy. This 
is useful in identifying galaxies with young stellar pop u- 
lations of ages up to 1 Gyr (see also Wil d et al.lkoiof ) 1 . 

In Figure [T] we show the distribution of the SDSS main 
galaxy sample in greyscale as well as contours in the 
PCI vs. PC2 plane. The locations of the starburst, 
post-starburst, passive, and typical star- forming galaxy 
populations are indicated. The colored dots traces out 
a t ime sequence for strong starbursts (details described 
by IWild et al.l l2010f ). These bursts can be fit with a 
model having an exponentially declining star formation 
rate with an e-folding time of 200-350 Myr and typical 
burst mass fractions of 10 to 25% of the stellar mass. To 
find a sample of recent /on-going starbursts we have then 
searched the SDSS DR7 Main galaxy sample for galaxies 
lying along or near this sequence. 

The second step in generating our sample was cross- 
correlating this starburst sample with the Galaxy Evo- 
lution Explorer (GALEX) far-ultraviolet (FUV) bright 
QSO catalog (GR3) to find QSO-starburst pairs. We 
used the following criteria: 

1. The background QSO is brighter than FUV = 19.5 

mag or F ° = 1.02 x 10 -15 erg cm -2 s _2 A _1 : 

& A=1300A to ' 

2. Signal-to- noise ratio > 8 in the SDSS g-band spec- 
trum of the galaxy; 

3. Galaxy stellar mass, M* > 1O 1O M 2 . We are in- 
terested in galaxies in the mass range where the 
galaxy population transitions from predominantly 
low- mass star-forming to high-ma ss quiescent sys- 
tems (e.g. iKauffmann et al.ll2003l ): 

1 PCA catalogs for SDSS DR7 are available for download at 
http : //www- star . st-and . ac . uk/$\sim$vw8/downloads/DR7PCA . html 

' z Stel lar mass catalogs, derive d from 5 band SDSS photom- 
etry b y Brinchmann et al. ( 2004) following t he prescription by 
Kauffmann et al. (2003) and Salim et al. (2007). The catalogs were 
obtained from http://www.mpa-garching.mpg.de/SDSS/ 
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Fig. 1. — Figure shows the first and the second principle co mponents, PCI an d PC2 respectively, from the Principal Component Analysis 
(PC A) study of SDSS spectra (rest-frame 3750-4150 A) by IWild etaLl ([2007). Contours show the distribution of SDSS galaxies in the 
PCI versus PC2 diagram and the greyscale shows their numbers in logarithmic units. The first principal component (PCI) provides the 
same information as the traditional D n (4000) index whereas the second principal component (PC2) provides additional information on any 
excess (or deficit) of Balmer absorption, over that expected for the 4000 A break strength of the galaxy. Therefore, a combination of the 
first two compo nents can id e ntify young stellar populations and hence starburst and post-starburst galaxies. The trajectory of starbursts 
from Figure 1 of Wild et al. (2010) is shown as colored dots indicative of their ages. The regions of the plot associated with various classes 
of galaxies with different star-formation histories are labeled in violet on the diagram. The galaxies from our sample are over-plotted as red 
symbols with their shapes representing their classification- diamonds for starburst and post-starburst, squares and circles for the control 
sample consisting of normal star-forming galaxies and passive galaxies respectively. 



Impact parameter to nearest UV bright QSO, p < 
200 kpc. This distance is approximately the virial 

1/3 



radius, R v i r = 



GM h 



215 kpc for galaxies 
10 12 M(7) (approximate 



_9tt 2 H^ 

with halo masses, Mh ~ iir~ iv±© 
halo mass for galaxy of stellar mass 10 10,5 M ); 

5. The foreground galaxy does not belong to a cluster 
as defined b y the SPSS galaxy cluster catalog (see 
iMiller et alJl2QQ5h , We are interested in the CGM 
in typical "field" galaxies. 

This resulted in six QSO-starburst pairs (hereafter 
the "starburst sample - SB/PSB"). We then repeated the 
above selection procedure for non-starbursts or "control" 
galaxies in SDSS DR7. Since there were many eligible 
QSO-galaxy pairs, we selected our targets so that we 
sampled the rest of the PCI vs. PC2 plane uniformly 
with galaxies that roughly has similar stellar masses and 
impact parameters as the starburst sample. For multiple 
pairs with similar properties, the one with the brightest 
background QSO in FUV was selected. This resulted in 
14 control galaxies with half of them passive (red) galax- 
ies with PC1> -1 and other half consisting of normal 
star- forming galaxies with PCI ranging between -3 and 
-1. The starburst and control samples are overplotted in 
the PCI versus PC2 plane in Figure [TJ while Figure [2] 
shows the SDSS multicolor image of the 20 target galax- 
ies from our sample. The target galaxy is at the center 
of the image and the direction of the QSO is shown with 
a blue arrow. The galaxy names are labeled in the top 
right corner and the ones where we have C IV detections 
are labeled in red. The symbols before the labels refer 
to the category of galaxies. We use diamonds to identify 



the SB/PSB sample and circles and squares for the con- 
trol sample of passive and normal star-forming galaxies 
respectively. 

The properties of the galaxies and QSOs in our sample 
are summarized in Table [TJ The name, position, red- 
shift of the foreground galaxy and the background QSOs 
are presented in columns 1, 2, 3, 4, 11, 12, 13, and 14 re- 
spectively. The impact parameter of the QSO sight line is 
presented in the last column. The PCI, PC2 and the cor- 
responding class of the foreground galaxy are presented 
in columns 5-7. Column 8 presents the stellar masses ob- 
tained from the MPA-JHU catalog, These masses were 
estima ted using methods described bv lBrinchmann et al.l 
(|2004f ). The orientations of the sightlines are presented in 
column 9, which represents the angle of the QSO sight- 
line with respect to the major axis of the galaxy. For 
example, any sight line along the major axis has O = 0°, 
whereas a sightline probing the halo perpendicular to the 
disk has O = 90°. The orientation was estimated using 
SDSS photometric i-band measurements. We considered 
all galaxies where the ratio of the major to minor axis 
(deVAB) > 0.7 to be face-on, and defined © as 0°. 

The star formation properties of the galaxies are pre- 
sented in Table [21 The Ha fluxes, luminosities, and ioniz- 
ing photons per second, Q* are presented in columns 2,3, 
and 4 respectively. These values were derive d by apply- 
ing d ust corrections (similar to that used by IWild et all 
2007, for emission lines) to the fluxes from the MPA-JHU 
catalog 3 and then converting to L(Ha) (see Eq.9.5 from 
iConti et al.l [2008). In this paper we wanted to evaluate 
the energy input to the CGM over the lifetime of the star- 

3 Line fluxes were obtained from the MPA-JHU catalog. 



4 



Borthakur et al. 



U075622 




1 J085252 



I J084356 





I J100801 



U101008 



• 



♦ J102846 



r 




♦ J122115 



♦ J122534 



i 



♦ J132107 




I J154527 





Fig. 2. — SDSS composite image of 37 .5 /; X for our sample of 20 foreground galaxies with stellar mass, M* between lO lo l -lO n - 4 M . 

At the median redshift of our sample of z = 0.09, the scale of the image corresponds to 63 kpc. The galaxies are identified by their name on 
the top right corner and the symbol represents their subclass: diamond for SB/PSB, square for normal star forming, and circle for passive 
galaxies. The red label indicate galaxies where C IV was detected. The blue arrow marks the direction of the background QSO. Details of 
our sample are presented in Table [T] 



burst (rather than using the instantaneous value). For 
this reason we do not measure the star formation rates 
for the starbursts using Ha. Instead, for the starburst 
sample we use the burst mass or burst mass fraction and 
the burst age to estimate a SFR or specific SFR (respec- 
tively) averaged over the over the duration of the burst 
(column 5). For t he control sample we use the SFR and 
specific SFR from iBrinchmann et al.l (|2004f ) (column 6). 
Both these estimates are for the central 3" region of the 



galaxies which is covered by the SDSS fiber. 

2.2. Cosmic Origin Spectrograph Observations 

The sample was observed with COS using segments 
B and A of the grating G140L, providing a resolution, 
R - 2000 - 3000 (a = 53 ± 10 km s" 1 ). The data were 
reduced using the standard COS pipeline and covered the 
wavelength range between rest-frame Lyman a to C IV 
with the exception of one control galaxy, J085254. Part 
of its spectrum (up to X re st < 1245 A) was missed as a 
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Fig. 3. — Rest frame Lyman a absorbers associated with 18/20 foreground galaxies for which we have data at the observed wavelength 
of the Lyman a . The dotted line shows the redshift of the target galaxy. The data are labeled with the target galaxy name in the upper 
right corner and its type is shown as a symbol next to the name: diamonds represent the 5 starburst galaxies, squares and circles represent 
the control sample of 6 normal star forming galaxies and 7 passive galaxies respectively. The best fit Voigt profile to the absorption feature 
is shown in red. 



result of an error in the setup of the grating. For the 
remaining 19 targets the data covered the entire range. 
The wavelength covered includes transitions of various 
other species such as C II, Si II, Si III, and Si IV. This 
allows us to probe both neutral gas via low-ionization 
transitions like Si II and C II as well as highly ionized 
gas via C IV. 

Although we have COS coverage for all the other sight- 
lines, we do miss a few transitions due to intervening 
absorbers or hot pixels at the position of the expected 
transition. For example, in the case of J140502, the 
saturated Lyman f3 transition of the background QSO 
completely absorbed the QSO flux at the expected po- 



sition of the Lyman a line associated with the fore- 
ground galaxy. Similarly, we have three galaxies namely 
J075622, J161913, and J132107, where the C IV is 
blended with the Milky Way's Al III, Lyman a of the 
background QSO, and lies at the position of noisy (hot) 
pixels respectively. No measurements could be made in 
these cases. This leaves us with a total of 18 targets with 
useful data at the expected wavelength of Lyman a and 
17 targets for C IV. 

We detected H I Lyman a in 15/18 targets. We also 
have unambiguously detection of C IV in 6/17 systems. 
Information on the properties of the Lyman a and C IV 
absorbers such as equivalent widths, full width at half 
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Fig. 4. — Rest frame C IV absorption features associated with our sample of 20 foreground galaxies. The data are labeled with the 
name of the foreground galaxy in the upper right corner and the galaxy-type is shown by the symbol. The fits to the C IV AA1548,1550 A 
doublet are plotted in red. The dotted lines show the position of the C IV doublet at the rest-frame of the target galaxy. The targets 
where no measurements could be made are indicated. 



max (FWHMs), and column densities are presented in 
Table 03 Figures [3] and H] show the Lyman a and C IV fea- 
tures at the rest frame of the target galaxy, also marked 
with a dotted line. The Lyman a lines are resolved and 
unsaturated barring a few cases. Two cases where the 
Lyman a lines are unresolved are marked in Table [3l 
The C IV doublets are also resolved and are unsatu- 
rated, based on the C IV A1548 to A1550 doublet ratio 
(which ranges from 1.54 to 1.87). The corresponding op- 
tical depths in the (weaker) C IVA1559 lines are ~ 0.2 
to 1. Apart from Lyman a and C IV, we do not detect 
any other major metal transitions except in the cases of 
J085301 and J092844. 



The absorption systems associated with target galaxies 
J085301 and J092844 showed multiple Lyman series and 
metal transitions. The absorbers are shown in Figures [5] 
& [6] and the measurements are presented in Table [U Be- 
sides multiple Lyman series transitions, these absorbing 
systems show the presence of low-ionization transitions 
such as Si II and C II as well as higher ionization species 
as traced by C IV. In J085301, the C IV A1548 transition 
lies at the edge of the COS spectral grating and the C IV 
A1550 transition falls in the gap. Hence in this case, the 
measurements presented in Table [4] were derived from 
the data with incomplete wavelength coverage and may 
also suffer from flux inaccuracies sometimes seen near the 
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Fig. 5. — Plots showing various transitions detected in the fore- 
ground galaxy J085301. Lines include H I, low-ionization metal 
transitions such as C II, and Si II, intermediate-ionization transi- 
tion like C III and Si III, along with hot gas traced by Si IV and 
C IV respectively. The Voigt fits to the profile are shown in red. 
We do not fit C III and Lyman 7 as they lie in the segment B of 
G140L grating where wavelength calibration is highly uncertain. 

edge of the gratings. 

We fitted Voigt profiles to the data to estimate all 
the parameters. The associated uncertainties were es- 
timated using the error analysis method published by 
ISembach fc Savage] (pj)92). Errors estimated using this 
procedure include continuum placement uncertainties, 
Poisson noise fluctuations, and zero point uncertainties. 
The continuum was established using absorption free re- 
gions of the spectrum within ±2000-4000 kms -1 of the 
absorption feature and subsequently fitted with a Leg- 
endre polynomial of o rder between 1 and 5 , similar to 
the procedure used by ISembach et all (|2004l ). This en- 
sures that the continuum is described in the neighbor- 
hood of the absorption feature precisely and nullifies any 
slow variations in the shape of the QSO's spectrum. The 
measurement also takes i nto account the approp riate line 
spread functions (LSFs; lOsterman et al.l (|201l[ )) for the 
aperture of the spectrog raph from the COS Instrument 
Handbook (|Dixonl 120101 ). For the non-detections, we 
present the limiting equivalent width, which represents 
the 3<j error in the equivalent width over 200 km s _1 i.e. 
+100 to —100 km s _1 in the rest frame of the foreground 
galaxy. Limiting column densities are estimated from the 
limiting equivalent widths using the relationship between 
the two for the linear part of the curve of growth. 
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Fig. 6. — Similar plots as in Figure but for the foreground 
galaxy J092844. 

We have corrected for inaccuracy in the absolute wave- 
length calibration of the COS data using Milky Way's 
strong ISM transitions. Since we can estimate the off- 
set of these transitions from the expected velocity with 
reasonable accuracy, we can then use them to correct 
for any offsets in the COS wavelength calibration. In 
order to find the intrinsic velocity offset of the ISM 
transitions, we use HI 21 cm measurements from 
the L eiden/ Argentine /Bon n (LAB) Survey of galactic 
H I (iKalberla et al.l [2QQ5h from the same region of the 
sky. Although the beam-size of the H I measurements is 
enormous as compared to the COS aperture, nevertheless 
we do not expect this would cause a large error as our ve- 
locity resolution is of the order ~ 100 km s _1 (full width 
at half maximum, FWHM). In most cases for the seg- 
ment A, we found velocity corrections to be of the order 
of 100 km s _1 . However, in the absence of ISM tran- 
sitions in the segment B, we could not apply the same 
corrections and the velocity calibration remains highly 
uncertain. In cases of multiple Lyman series transitions 
located in both the segments, the ones in segment B (e.g. 
Lyman 7) show a velocity offset of ~ —700 km s _1 rela- 
tive to those in segment A (Lyman a ). 

3. RESULTS 

3.1. The Strength of Lyman a Absorbers 

The properties of the CGM as a function of properties 
of the host galaxy are not yet clearly understood. Studies 
have detected gas traced by Lyman a out to radii of a 
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Fig. 7. — Plot showing the variation of the Lyman a equivalent widths as a function of the impact parameter from the foreground 
galaxies. Lyman a data for our sample are sho wn in red with sym bols i ndicating galaxy class ( same as in Figure^). The data in blue and 
green are fro m the recen t ly pub lished work bv lThom et al.l (120121 ) and [Prochaska et al. (2011J) respectively. The magenta line shows the 
best fit from lChen et al.l (^OOlbT ) to their Lyman a absorber sample. 
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Lyman a equivalent width is plotted against the ori- 
entation of the sight lines, defined as the angle between a sightline 
and the major axis of the target galaxy. The orientation is 0° when 
the sightline probes the CGM along the major axis and 90° when 
the sightline is along the minor axis. All data points with an orien- 
tation of exactly 0° are for face-on systems. The symbols indicate 
galaxy types as in previous figures. 



few h undred kpc dChen et al.l ll998l l2001ri iBowen et all 
120021 ; ISteidel et al.l 120 111 and other authors) and found 
a weak inverse correlation between the strength of the 
Lyman a line and the impact parameter. 

Figure [71 shows Lyman a equivalent width versus im- 
pact parameter. Our data are shown as red symbols with 
their shapes representing their spectral classes based on 
the indices PCI and PC 2. The m a genta l ine shows 
the best fit obtained by IChen et all ([2001bf ) for their 
data. We also overplot data from recent studies by 



iProchaska et al.l ([201 lh and lThom et al.l (|2012l) for a vi- 
sual comparison. We do not see a trend of equivalent 
width with impact parameter in our sample. This is not 
surprising since our sample spans a rather narrow range 
in impact parameter. Furthermore, these objects cover 
the full range of equivalent width seen in other surveys. 
Interestingly, the top five strongest absorbers are associ- 
ated with the CGM of the control sample. 

Some authors have found a correlation between the 
strength of Lym an a absorption and the luminosity o f 
the host galaxy ([Chen et al.ll2001bl : IBowen et ani2002l ). 
We see no correlation between the Lyman a equivalent 
width and the galaxy stellar mass, which is a proxy for its 
luminosity. This may be because we span a rather narrow 
range in stellar mass (10 101 — 10 11,4 M ), although the 
implied range in the dark halo mass is considerably larger 
(few xlO 11 to a few xlO 13 M ). 

We also investigate the strength of the Lyman a tran- 
sition as a function of orientation of the sightline with re- 
spect to the galaxy major axis (Figure [8]). The strongest 
absorbers are aligned close to the major axis, although 
some of the weak absorbers can also be found there. 
However, due to small number statistics of our data a 
correlation cannot be verified with certainty (the Spear- 
man's rank analysis implies that there is a ^20% chance 
that there is no correlation). 

A more intuitive way of understanding the com- 
bined effect of orientation and impact parameter on the 
strength of the absorbers is presented in Figure [U The 
figure shows a cartoon of a generic galaxy at the center. 
The equidistant regions from the center are marked with 
dotted semi-circles. It is to be noted that this model is 
just for illustration and does not represent the range of 
morphologies and sizes of our target galaxies. In order to 
distinguish between the two sub-samples, the starburst 
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Fig. 9. — A graphical representation showing our sightlines with respect to a model galaxy. The plot shows impact parameters, 
orientations, and galaxy-types, as well as the strength of associated Lyman a and C IV absorbers for each galaxy in our sample. All the 
data points with an orientation of exactly 0° are for face-on systems. For easy comparison the starburst galaxies are plotted on the right 
side of the plot and the control sample consisting of the normal and passive galaxies on the left. The dotted circles mark the loci of impact 
parameter of 100, 150, and 200 kpc from the center of the galaxy. The strength of the absorber is coded as the size of the symbol. Note 
that the symbol size scale is not linear. 
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Fig. 10. — Plots showing the variation of Lyman a equivalent width as a function of PCI, which represents the 4000 A break (D n (4000)) 
and as a function of the specific star-formation rate (SFR/M* in inverse years). The galaxy type is represented by the symbols as in 
previous figures. 



galaxies are plotted on the right side of the figure and 
the control galaxies containing the normal and passive 
galaxies are plotted on the left. It is evident that the 
Lyman a absorbers with equivalent width WL ya > 1A 
are only found in the control sample. 

To further emphasize the difference in the average Ly- 
man a strength between the starburst sample and the 
control sample, we study the variation in terms of both 
the value of PCI and the star- formation rate per unit 
stellar mass (the specific star- formation rate see Table [2]). 
Figure [T0l suggests that the scatter in the Lyman a equiv- 



alent width is larger in the passive galaxies and smallest 
in the starbursts. 

3.2. Highly Ionized Gas in the CGM in Starburst 
Galaxies 

A drop in the strength of the Lyman a line can be due 
to either a relative lack of gas in the CGM of starburst 
galaxies or to a higher degree of ionization and/or tem- 
perature. A distinct difference between the two scenarios 
would be an enhanced presence of absorbers associated 
with highly ionized and/or hot (coronal-phase) gas. The 
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Fig. 11. — Histogram showing the rate of C IV detection (and 
non detections) as a function of PCI and the specific SFR respec- 
tively. The SB/PSB and the control samples have been labeled. 
These figures show the higher detection fraction for SB/PSB sam- 
ple as compared to the control sample. Chi-squared contingency 
test confirms that the two samples have different C IV detection 
rates with 99.9% confidence. 

best probe of such gas in our data is the C IV line. Our 
COS data provided useful information on this line in 17 
cases (5 starbursts and 12 control galaxies). It is then 
possible to look for any correlation between the presence 
of a starburst and of highly ionized gas traced by C IV. 

We detected C IV absorption- lines in six galaxies. 
As described above, we fit Voigt profiles to the dou- 
blet and this yielded typical column densities of several 
xlO 14 cm -2 (see Table [3]). Four of the detected galaxies 
are starbursts and the remaining two are control galax- 
ies. Thus, we find C IV absorbers in 80% (4/5) of the 
SB/PSB and in only 17% (2/12) of the control galaxies. 
Interestingly, the two detections of the C IV absorbers in 
the control sample are associated with passive galaxies 
and none with normal star forming galaxies (see Fig- 
ured!]). 

Figure [12] shows the equivalent width as a function of 
impact parameter and specific SFR. By design, we do 
not explore a large impact parameter range and hence 



no tren ds are prominent unlike that found bv lChen et al.l 
(|2001al ) for a sample of galaxies at a median redshift of 
0.4. Our da ta (plotted in red) is combined with that of 
I Chen et al.l (plotted in blue) is shown in the left panel 
of the figur e . Follo wing the fitted empirical relation in 
IChen et all (|2001a[ ). we have scaled the impact param- 
eters for our sample by the factor (M/M*) -0,51 . Here 
M is the total galaxy stellar mass and M* is the fidu- 
cial value for a Schechter function fit to the galaxy stel- 
lar mass function: M* = 7 x 10 10 M (e.g. iLi fc White! 
l2QQ9h . This panel clearly reinforces our previous con- 
clusion that the CGM surrounding the starbursts has 
unusually strong C IV absorption-lines at large impact 
parameters compared to typical galaxies. We do ex- 
plore a large parameter space in terms of specific SFR 
and find strong absorbers in galaxies with high specific 
SFRs. A qualitatively sim ilar trend was also observed by 
iTumlinson et al.l (|2011b[ ) for their O VI absorber sample 
(plotted in blue in the right panel of figure [12]) . 

Recent numerical simulations have modeled CGM 
properties in the presence of galact ic outflows using 
SPH-based cosmological simulations (Fo rd et "all l2Q12h 
and adaptive-mesh refi nement simulat i ons f or Milky- 
Way mass e d gal axies ([Hummels et al.l l2012[ ). While 
iFord et all (|2012f ) found a subset of C IV absorbers at 
impact parameters of 100 kpc from their host galaxy 
that had column densities greater than 10 14 cm -2 , the 
median value was about ~ 10 13 cm -2 (corresponding 
to a C IV equivalent width of 0.05 A). Simulatio ns of 
Milky- Way- like galaxies by iHummels et al.l (|2012f ) pre- 
dicted typical C IV columns at impact parameters > 100 
kpc of r^j 10 12 cm -2 . Both sets of simulations are consis- 
tent with our data for typical galaxies, but fall far short of 
the observed C IV column densities of a few x 10 14 cm -2 
in the CGM of starbursts. 

We note that most of the spectra in our sample have 
similar signal-to-noise ratios and therefore the C IV non- 
detections are not a consequence of lower data quality in 
these spectra. Figure [13] shows the stack of all the non- 
detections in the two control sub-samples. This was done 
in order to detect any weak signature of the presence of 
C IV absorption in the data by increasing the S/N of the 
averaged spectrum. The C IV spectrum of each galaxy is 
stacked at the peak velocity of associated Lyman a after 
correcting for the intrinsic COS wavelength calibration 
corrections (as described in Section 12^2]) . The wavelength 
corrections help against dilution of the stacked spectra. 
The resulting uncertainties in rest-frame velocities are 
< 100 km s _1 , which is lower than the resolution of the 
grating or the FWHM of all the C IV detections. The 
stacked spectra were generated by adding all the C IV 
non-detections for the two control galaxy subsamples i.e. 
5 normal galaxies and 5 passive galaxies each. The single 
case of C IV non-detection in the starburst sample is not 
shown here. The limiting equivalent width of C IV for 
the normal star- forming galaxy stack is < 118.5 mA (3<r), 
compared to values of ~ 400 to 900 mA in the six 
detections. Passive galaxies do show a slight dip, however 
due to undulating nature of the noise in the spectrum this 
feature is still within the noise. So, in neither stacked 
spectrum is there a significant absorption feature. 

There are important differences between C IV ab- 
sorbers in starburst and control galaxies besides the de- 
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Fig. 12. — Plots sho wing the st r ength of C IV as a function of normalized impact parameter (left) and specific SFR (right). The 
scaling (X/X*) used by Chen et al. ( 200 1a) was in terms of luminosity i.e. L/L*. Since our ta rgets are starburst galaxies with very high 
luminosities, we choose to adopt a stellar mass scaling which is more appropriate to Chen et al.'s scaling. Here, M is the tota l galax y stellar 
mass and M* = 7 x 10 10 Mq. The three galaxy classes are l abeled. In green we overplot data on C IV from lChen et al.l (|2001aD on the 
right and in blue data on O VI from Tumlinson et al. (2011b) on the left. 

These spectra were obtained by stacking both detections 
and non-detections together for each of the transitions - 
Lyman a A1216 A, C II A1334 A, C IVAA1548,1550 A, 
Si III A1206 A, and Si IV AA1393,1402 A. The dotted 
lines show the rest-frame of the absorbers based on the 
position of the Lyman a centroid. The Lyman a spectra 
were obtained by adding each individual spectrum cor- 
rected to the rest-frame velocity of the target galaxy. For 
the other transitions that are weaker than Lyman a, we 
correct the velocity to the rest-frame Lyman a velocity 
whenever Lyman a was detected. If not, we used the 
optical redshift of the target galaxy to set the rest-frame 
velocity. 

Lyman a was clearly detected in the stacked compos- 
ite spectra in all sub-samples and C IV was detected for 
starburst and control-passive galaxies. We also detected 
Si III in all three composite spectra at a strength of about 
20-25% of Lyman a. Table [5] summarizes the equivalent 
widths, FWHMs, and column densities of all the detected 
transitions and provides upper limits for non-detections. 
Among the detections, the Lyman a equivalent width as 
well as its velocity spread increased from the starburst to 
normal to passive galaxies. It is not surprising that the 
equivalent width and the velocity spread of the compos- 
ite Lyman a features are correlated. This is due to the 
contribution of strong (saturated) Lyman a lines associ- 
ated with some of the passive and normal galaxies. Since 
the stacked spectra were produced by averaging both sat- 
urated lines with unsaturated and even non-detections, 
the resulting composite may not look saturated. 

Another way to study the difference in the ionization 
state of the CGM in the starburst halos is shown in Fig- 
ure [15l Here we plot the ratio of the equivalent widths 
of C IV to Lyman a as a function of PCI and specific 
star formation rate. The top panels show the ratios for 
individual galaxies in our sample. Not all galaxies have 
detectable Lyman a and the plot includes only those with 
Lyman a detections. The upper limits in the plot rep- 
resent cases where no C IV was detected. None of the 
four galaxies with an upper limit on the Lyman a equiv- 
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Fig. 13. — Composite rest-frame C IV spectra obtained by stack- 
ing all non-detections for normal star-forming and passive galaxies. 
The stacks contain 5 galaxies each. The dotted lines mark the ex- 
pected position of the C IV doublets. The A1548 A feature is at 
zero velocity and the A1550 A feature is at 498 km s _1 . 

tection rate. While the presence of C IV in the con- 
trol (passive) galaxies is accompanied by the presence of 
strong low-ionization metal transitions, the same is not 
true for the starburst galaxies. The two passive galaxies 
with C IV detections show transitions such as Si II, C II, 
Si III, C III, and Si IV. The equivalent widths of the Ly- 
man a absorbers in each of these two systems (1.5 and 
2.9 A) are also significantly higher than those found in 
the starburst galaxies (0.3 to 1.1 A; median of 0.5 A). 

To investigate this more generally, we have created 
stacked composite spectra for the starburst and control 
samples. The composite spectra for Lyman a , C II, 
C IV, Si III, and Si IV for the SB/PSB sample and 
the two control sub-samples are presented in Figure [Ml 
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Fig. 14. — Composite rest-frame spectra of galaxies from the three classes. Each composite was obtained by stacking all the non-corrupted 
spectra including detections and non-detections. The panels show Lyman a , C IIA1334, the C IV doublet, Si IIIA1206, and Si IVA1393. 
The dotted line(s) marks the expected rest-frame position of the transitions. The Voigt profile fits to the spectra are shown in red solid 
lines. For non-detections, we plot red dashed lines showing the limiting equivalent width. 



alent widths show C IV absorption. In other words, we 
have detected Lyman a in all the galaxies where we have 
detected C IV, with the possible exception of J140502, 
where we have no information on the Lyman a transition. 

The bottom panels of Figure [15] show the subset of 
galaxies where we have detected C IV and also include 
the average C IV/Lyman a ratio for the three sub- 
samples in green. The averages were derived from the 
composite spectra obtained by stacking the velocity cor- 
rected spectrum of individual galaxies. The larger mean 
ratio of C IV/Lyman a in the starbursts is evident. This 
would suggest that the CGM has a higher overall degree 
of ionization in the starbursts. We believe that this is 
an ionization effect, and not due to higher metallicities 
in the CGM of the starburst galaxies. If it were just a 
metallicity effect, then the starburst galaxies should show 
strong absorption-lines from metals in the lower ioniza- 



tion states in addition to C IV. This is not the case as 
can be seen from our data (see Table [5]). 

4. DISCUSSION 

The results presented above (section 3.2 and Ta- 
bles [3] and [5]) show that material with a C IV column 
density of a few xlO 14 cm -2 covers a large fraction (~ 
80%) of the sight-lines through the halos of starburst 
galaxies extending out to radii of roughly 200 kpc. 

To get a crude estimate of the mass represented b y the 
C IV absorbers we follow [Tumlinson et al.l (|2011b| ) and 
use a simple model-independent geometric method. We 
take the total number of C IV ions in the CGM to be 
given by the product of the C IV column density for our 
absorbers, the cross-sectional area of the CGM that we 
probe (7rr 2 , where r = 200 kpc), and the covering fac- 
tor of the C IV absorbers within this radius (4/5). This 
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Fig. 15. — Plots showing the ratio of C IV to Lyman a equivalent widths as a function of PCI in the top left panel and specific SFR in 
the top right panel. These plots include galaxies for which we have associated Lyman a absorbers. That is, galaxies with upper limits to 
their Lyman a equivalent widths or those where we do not have measurements of either Lyman a or C IV are not included in these plots. 
In the bottom p ane l, we omit the individual upper limits and instead plot the average for each of the classes using the composite spectra 
shown in Figure fl4l 



implies a total mass of C IV within a radius of 200 kpc 
of roughly a few xlO 6 M . For collisionally ionized gas 
(see below), the minimum ionization correction (which 
occurs for a temperature of 10 5 K) implies a total Car - 
bon mass of about 10 7 M (jDopita Sut her land fl996h . 
The maximum likely metallicity of the CGM would be 
solar. In this case the minimum total gas mass associ- 
ated with the C IV absorbers would be ~ 3 x 10 9 M 
(but is likely to be larger). Since the mean stellar mass 
of the starbursts in our sample is about 2 x 10 10 M , 
we conclude that the C IV absorbers trace a significant 
repository of the baryons in these galaxies. 

The results presented above also suggest that star- 
bursts are affecting the ionization state of a significant 
amount of gas in the circum-galactic medium. We will 
first discuss whether this is physically plausible, will then 
use these results to estimate the properties of the absorb- 
ing material, and then briefly discuss the implications of 
these results. 



4.1. Photoionization 



One possibility is that the C IV in the starburst galax- 
ies arises in photoionized gas. The radiation field in- 
cident on a CGM cloud of density, n, at a distance, r, 
from the galaxy can be characterized using the ionization 
parameter, U . that includes contribut ion from both the 
extragalactic Haardt fc Madaul ([2012) background and 
the ionizing flux from the starburst: 



U = Uhm + U SB 



(1) 

As described by iTumlinson et all ([2011a[h U can be ex- 
pressed as 

- //js , Q*ft 

\®HM 



u 



c n 



47rr 2 



(2) 



where Q* is the total rate of hydrogen ionizing photons 
produced by the starburst with an escaping fraction of 
/esc, t is the distance from the starburst to the cloud 
(which we approximate as the impact parameter) and c 
represents the speed of light. 

We present two scenarios for photoioniza tion of a CGM 
cloud as modeled using the code Cloudy (|Fer land et al.l 
Il998f ) for a plane-parallel slab. These two cases can be 
considered to bound the range of possibilities. These sim- 
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Fig. 16. — Photoionization models using Cloudy for a CGM cloud that is ionized by the diffuse metagalactic ionizing background. 
Different species are labeled using different colors. For easy viewing the C IV and O VI transitions are additionally marked with filled 
circles and filled squares respectively. The thick dashed line represents the H I column density and was set as a termination criterion for 
the simulations. The right panel shows the ratio of the various species in terms of C IV. These ratio s are i ndependent of the termination 
criterion. The metagalactic flux at redshift z=0 is normalized to that described bv lHaardt fc M adau (20J3) • 

bust matrix for estimating the ionization parameter. The 
plot in the right panel of Figure [T6l shows these ratios nor- 
malized to C IV. The strongest constraint on the ioniza- 
tion parameter comes from the lower limit on the ratio of 
the Si IV and C IV columns. This requires log(U) > -1.5. 
We use a value for $hm = 2750 photons cm -2 s _1 based 
on the most recent mod eling of the UV background by 
lHaardt fe Madaul (120121). consist ent with the value found 
earlier bv iGiroux fe ShulU (|1997[ ). We then increase this 
flux by a factor of two to account for illumination of both 
sides of the photoionized slab. This yields an upper limit 
to the density of the cloud n < 6 x 10 -6 cm ~ 3 . We find 
a minimum ionization correction to go from C IV to C is 
a factor of four. To produce the observed C IV column 
density would then require a cloud size > 240 [Z0/Z c i ou d] 
kpc (where the scaling accounts for the the cloud metal- 
licity). For plausible (sub-solar) metallicity the implied 
sizes are significantly larger than the CGM. Moreover, 
photoionization by the metagalactic background would 
not naturally explain why strong C IV is much more 
common in the CGM of starbursts than in the control 
sample. We conclude that this model is implausible. 

Next, we explore photoionization of the CGM caused 
by the radiation escaping from the starburst galaxies. 
Figure [TT] shows the distribution of various species in a 
CGM cloud in terms of the ionization parameter pro- 
duced by a starburst. The starburst spectrum used for 
these runs is a model spectrum for a starburst with 
constant sta r formation for 100 M yrs as produced by 
Starburst 99 (|Leitherer et al.l 12010). Because the ioniz- 
ing spectrum of a starburst is softer than that of the 
metagalactic background, the constraints on the ioniza- 
tion parameter are even more severe. The lower limit to 
the ratio of Si IV and C IV columns implies log(U) > 
0.2. We estimate Q* using the extinction corrected Ha 
fluxes and the values are presented in Table [2j Adopt- 
ing the mean values for (4 x 10 53 s _1 ) and r (160 
kpc) for the starburst sample, the required density is 
n < 2 x 10 -6 f esc cm -3 and the implied cloud size is 
> 130 [Nc/Nciv][Z /Z c ioud]felc k P c - Observations of 



Fig. 17. — Photoionization models using Cloudy for a CGM cloud 
that is ionized by the starburst in the host galaxy. Different ionic 
species are color coded. The C IV and O VI transitions are addi- 
tionally marked with filled circles and filled squares, respectively. 
The spectrum of the starburst is modeled using Starburst99 for a 
100 Myrs old continuously star forming galaxy. The plot shows the 
ratio of column densities of various species with respect to column 
density of C IV as a function of the ionization parameter, Usb> 
which can be expressed in terms of the physical parameters of the 
cloud and the starburst as described by Eq. [2] 



ulations assume the relative abundances of the different 
metals have solar values with no dust depletion. 
In the first case, we use the Haardt fc Ma dau meta- 

1 120121 ). We illu- 



galactic background (jHaardt fc Madau 
minate the slab from two sides in this case to account 
for the isotropic nature of the metagalactic background. 
Figure [16] shows the ionization of the various species as 
a function of ionization parameter, which in turn is re- 
lated to the density of the cloud (see EqEJ. The left panel 
shows models integrated until an H I column density of 



10 



14.5 



cm is reached (matching our inferred values in 



the CGM of the starbursts). Changing this limit moves 
the plotted lines up and down such that their ratios re- 
main the same. Hence, their ratios can be used as a ro- 
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typical local starbursts yield an u pper limit on the esc ape 
fraction of f esc =0.01 or 1% (e.g. iGrimes et al.ll2QQ9h . so 
this scenario would yield cloud sizes even larger than in 
the case of the metagalactic background. 

4.2. Collisional Ionization & Galactic Winds 

If the C IV does not trace photoionized gas, the al- 
ternative is that the gas is hot and collisionally ionized. 
This process can not be due to shocks in the CGM as- 
sociated with a major galactic merger. The images of 
our starburst sample in Figure 2 show that they are 
mostly normal single gal axies. This is consistent with 
iSanders fc Mir abell ([ 19961 ). who find that only 12% of low 
luminosity starbursts similar to our targets are merger 
driven. Since our C IV detection rate ~80%, we there- 
fore conclude that any collisional ionization of the CGM 
is not primarily merger-driven. 

Let us then explore a model in which the C IV traces 
gas clouds in the CGM that are shock-ionized by a 
starburst-driven wind that may or may not be escaping 
(galactic fountains). We first consider whether such out- 
flows driven by starbursts with the ages and strengths 
corresponding to our sample would be able to affect 
th e CGM at s uch la rge radii. The models described 
in IWild et all (|201Q[ ) applied to our sample yield av- 
erage burst ages of 200 Myr, average burst masses of 
~ 10 9 M©, average burst mass fractions of 11%, average 
star-formation rates during the burst of ~ 7M yr _1 . 

The best-studied starburst wind is M82. O bserva- 
tions with Chandra ([Strickland fe Heckmanll2009[ ) imply 
that the outflow speed of the wind fluid produced in- 
side M82 is about 2000 km s _1 and the mass and kinetic 
energy outflow rates in this material are ~ 5 M yr _1 
and 7 x 10 43 erg s _1 respectively. Our estimated burst 
parameters imply an average star formation rate about 
70% as high in our starbursts as in M82 so we will 
scale the M82 outflow rates down by a factor of two for 
their winds. These empirical estimates of the wind out- 
flow rates are very similar to theoretical estimates for a 
starburst wi th a star-formation ra te of 7 M yr _1 (Star- 
burst 99; see iLeitherer et al.ll2010l ). 

The simplest model to consider is if the wind prop- 
agates freely outward at 2000 km/sec into the halo. 
The wind will then travel ~200 kpc in 10 8 years (con- 
sistent with our results). More realistically, we can 
consider a wind that propagates into a spherically- 
symmetric volume-filling CGM, thereby creating an ex- 
panding wind-blown bubble . Th is ge neral problem 
has be en analyzed by iDvsonl (|1989[ ) and iKoo fc McKed 
(fl992h for power-law radial density profiles. Taking an 
initial radial profile for the CGM corresponding to an 
isothermal halo (p oc r -2 ), the radius of the expanding 

• 1/3 1/3 

bubble is given by r = 60E 4 % M C J M 10 t$ kpc. Here E is 
the rate of kinetic energy injected by the starburst given 
in units of 10 43 erg s _1 , the CGM mass out to a radius 
of 200 kpc is in units of 10 10 M and the time since the 
starburst began is given in units of 10 8 years. 

The gas mass of the CGM is very uncertain, but is 
likely to be comparable to or a bit less than the galaxy 
stellar mass (e.g. Hummels et al. 2013). Since the av- 
erage stellar mass of the starbursts is 2 xlO 10 M , the 
fiducial value for Mcgm, io = 1 is reasonable. We take 
an average age of ~ 200 Myr and an average kinetic en- 



ergy injection rate of 5xl0 42 erg sec -1 for our starburst 
sample. This leads to an outer radius for the expanding 
wind-blown cavity in the halo of about 100 kpc. Given 
the highly idealized nature of our model and the uncer- 
tain mass of the CGM, we regard this simple estimate 
as only showing that it is plausible for a starburst-driven 
wind to affect the bulk of the CGM. 

Based on the properties of our Milky Way and of QALS 
associated with other galaxies, we know that the real 
CGM contains clouds or clumps. Assuming a simple pic- 
ture in which the CGM cont ains both clouds and a tenu- 
ous inter-cloud medium (e.g. iKaufmann et al.l l2009). the 
outflowing wind will then not only inflate a large bub- 
ble/cavity in the tenuous component, it will also shock- 
heat and accelerate the population of imbedded ambient 
CGM clouds. 

To produce substantial amounts of C IV in these 
clouds requires a post-shock temperature above 10 5 K 
and hence shock spee ds in excess of 100 km s _1 (e.g. 
IShull fc McKeel 119791 ). For a cloud shocked by a wind, 
the velocity of the cloud shock is related to the velocity 
of the wind by v shock = v wind (p wind / pdoud) 1 /2 , where 
pwind and pdoud are the density in the wind and the 
pre-shocked cloud. Our adopted wind parameters give 
Vwind = 2000 km s _1 and p wind = 4 x 10" 31 g cm -3 
at r = 150 kpc. This implies that a cloud density of 
Pdoud < 1.8 x 10 -28 g cm- 3 (n c i oud < 1.1 x lO^cm -3 ) 
would lead to v s h oc k > 100 km s _1 . 

Adopting the minimum corrections for ionization and 
metallicity discussed above, the observed C IV column 
densities imply minimum total hydrogen column densi- 
ties of at least a few xlO 18 cm -2 . Together with the 
cloud density above, the implied minimum cloud size 
would be ~ 10 kpc. This is much smaller than the total 
path- length our line-of-sight traverses though the CGM 
(roughly given by the impact parameter). This size mis- 
match is consistent with the absorbing material having 
a small volume filling factor (i.e. being clumpy). The 
associated time for the shock to traverse the absorbing 
cloud (the cloud crossing time) would be greater than 10 8 
years. Wind-shocked clouds can survive several cloud- 
cross ing times before be ing; hydro dyamically disrupted 
(e.g. iPittard et al.ll2010l h consistent with the continued 
existence of the shocked clouds over a period of a few 
hundred Myr (the starburst duration). 

These estimates can be compared to predictions for col- 
umn densities behind shocks. The total column density 
of material cooling behind a radiative shock is indepen- 
dent o f the gas density and is g iven by N rad = 10 17 5 v 4 00 
cm -2 ([Draine fc M cKee 1993). Detailed shock models 
yield predicted C IV col umn densities of ~ 1 13 cm~ 2 
for v shock ~ 100km s" 1 (jShull fc McKeel fl979h . Since 
the observed column densities are about an order-of- 
magnitude larger than this, it seems likely that the C IV 
absorption we see arises from an ensemble of multiple 
shocks. 

The observed breadth of the C IV lines (FWHM ~ 
120 to 320 km s _1 ) could be associated with wind-driven 
shocks. For a CGM cloud with a total gas column den- 
sity of 10 19 cm -2 located initially at rest at a radius 
of 150 kpc from the starburst, and then exposed to 
the ram pressure of our adopted wind, Equation 2 in 
iHeckman et al.l (|201l[ ) predicts the cloud terminal veloc- 
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ity would be about 200 km sec . 

Finally, we can consider the energetics of such a pic- 
ture. If the C IV traces gas that is cooling radiatively, 
this cooling rate should not exceed the total heating rate 
available from the starburst wind. For the masses, den- 
sities, and temperatures discussed above the total CGM 
cooling rate traced by the C IV absorbers wo uld be 
- 3 x 10 42 erg s" 1 ([Sutherland fc Dopital fl993l ) . This 
rate is independent of metallicity (a lower metallicity in- 
creases the total gas mass and thermal energy implied 
by the observed C IV column densities, but increases the 
radiative cooling time of the shocked gas by roughly the 
same amount). It is still a minimum rate since it as- 
sumes the gas is cooling from a temperature of 10 5 K 
(the minimum needed to produce C IV). This cooling 
rate is similar to the kinetic energy outflow rate in the 
wind that we have estimated above based on scaling the 
M82 wind (5 xlO 42 erg s _1 ). It is possible then to pro- 
vide the necessary heating, but a significant fraction of 
the wind kinetic energy is needed. 

4.3. Relationship to VI Absorbers 

In an important paper, iTumlinson et al.l (|2011bf ) re- 
ported on results obtained with COS that probed the 
O VI doublet (1032,1038 A) along sightlines through the 
CGM of 42 galaxies at redshifts of 0.1 to 0.36 and with 
impact parameters ranging up to 150 kpc. They find a 
very high detection rate (90%) of O VI in the normal star- 
forming galaxies (specific SFR ~ 10 -11 to 10 -9 ' 5 yr _1 ) 
and a much lower detection rate (33%) among the passive 
galaxies (specific SFR < 2 x 10 -12 yr _1 ). At impact pa- 
rameters similar to those we have probed (p > 100 kpc) 
the typical O VI column densities are ~ 10 14,5 cm -2 in 
the CGM of the detected galaxies. 

While both their results and ours link the presence of 
highly- ionized gas in the CGM to star formation in the 
central galaxy, there is an intriguing difference: we find 
that strong C IV absorbers are preferentially detected in 
the CGM of starbursts, but are not detected in normal 
star- forming galaxies. The weakness of the C IV absorp- 
tion at impact parameters like these is also seen in the 
Chen et al. (2001a) sample of normal disk galaxies at 
< z >~ 0.4. 

Putting all these results together would imply that the 
ratio of O VI and C IV column densities in the CGM of 
normal star forming galaxies is greater than ~ 10. This 
could imply the O VI in these galaxies traces exceed- 
ingly low density material that is photoionized by the 
metagalactic background. However, the Cloudy models 
discussed above would require densities less than 6 x 10 _7 
cm -3 in this case. The observed O VI column densities 
then lead to sizes greater than a Mpc (even for solar 
Oxygen abundances). This is not tenable for the CGM. 

Alternatively, the large ratio of O VI to C IV could in- 
dicate collisionally ionized gas significantly hotter than 
T = 10 5 K. We have argued that the C IV in the starburst 
CGM arises in collisionally ionized g as as well. Our re- 
sults that C IV is much stronger in the CGM of starbursts 
might suggest that the O VI lines would also be stronger 
in the CGM of starbursts than in normal star forming 
galaxies. Currently, the data to establish that result is 
lacking. Observations of O VI along site lines directly 
into starbursts similar to those in o ur sample show hig h 
column densities of ~ 10 15 cm -2 ([Grimes et al.ll2009| ). 



However it is not clear if this gas is actually located in the 
CGM (rather than close to the starburst). Future obser- 
vations of O VI in the CGM of starbursts with HST/COS 
would be most instructive. 

5. SUMMARY & IMPLICATIONS 

The circum-galactic medium (CGM) comprises of gas 
located in the halo of a galaxy extending out to the virial 
radius. As such, it is the interface between a galaxy and 
the inter-galactic medium (IGM). The effect of winds 
driven by intense star formation on gas in the CGM may 
be crucial to the evolution of both galaxies and the IGM. 
For example, can such negative feedback shut down the 
delivery of fresh gas flowing through the CGM into a 
galaxy and thereby quench future star formation? If so, 
does this play a role in the on-going migration of galaxies 
from the blue star-forming population to the red passive 
one? Can metals and energy be transported through 
the CGM and ultimately be delivered to the IGM? An- 
swering these questions requires detailed observations of 
galactic winds. The best-studied such winds are driven 
from starbursts in the nearby universe. While existing 
data demonstrate that such starburst-driven winds ex- 
ist, they provide scant information on their properties 
beyond radii of a few tens of kpc and so do not probe 
the bulk of the CGM. 

In order to probe the effect of starburst-driven winds 
on the CGM we conducted a COS UV spectroscopic 
study using background QSOs to probe the CGM (in ab- 
sorption) around foreground galaxies. This study probed 
sight-lines out to radii of 200 kpc for 20 galaxies. We tar- 
geted starburst galaxies and control galaxies consisting 
of normal star forming and passive galaxies. We com- 
pared the properties of CGM in these three classes and 
found the following: 

1. Most (83%) of the target galaxies in our sample 
showed Lyman a absorbers in the CGM. The disper- 
sion in Lyman a equivalent width was larger for passive 
galaxies, with some absorbers having values of 3 A and 
others with upper limits whose values were an order-of- 
magnitude lower. This is unlike the starburst sample 
where Lyman a was detected in all cases but had lower 
strengths (0.4 to 1.1 A). 

2. Highly ionized gas traced by C IV, was predom- 
inantly found in the CGM of starburst galaxies. The 
C IV line was detected in 80% of the starburst galaxies 
(4/5) compared to only 17% (2/12) of the control galax- 
ies. These lines are resolved (intrinsic FWMH ~ 100 to 
320 km s _1 and are mildly saturated with implied C IV 
column densities of a few xlO 14 cm -2 . 

3. Multiple transitions from low- ionization species 
such as Si II and C II as well as intermediate-ionization 
transitions like Si III and C III were detected in the two 
passive galaxies where C IV was detected. This was not 
the case for the four starbursts with detected C IV. This 
implies that the CGM in the starbursts is more highly 
ionized than the CGM in the passive galaxies, even those 
where C IV is detected. This inference is further sup- 
ported by the result that the ratio of the equivalent 
widths of the C IV and Lyman a lines is larger in the 
starburst galaxies than in either the passive or normal 
star forming galaxies. 

4. It is very unlikely that the C IV absorbers in the 
starbursts are photoionized by either the metagalactic 
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background or the stellar radiation from the starburst. 
CLOUDY models show that high ionization state of the 
gas requires cloud densities that are very small. To ex- 
plain the observed column densities requires path-lengths 
through the absorbing gas that are significantly larger 
than the size of the CGM. 

5. We have shown that a starburst with properties like 
those we infer for our sample members (average burst 
ages of ~200 Myr and average burst mass fractions of 
~10% of the stellar mass) can plausibly drive an outflow 
that can affect the bulk of the CGM. In such a model 
the C IV would arise in pre-existing clumpy clouds or 
filaments in the CGM that are being shock- ionized and 
accelerated by the ram pressure of the wind. The ensem- 
ble of such clouds must have a large areal covering factor 
(~ 80%). 

6. Even for the most conservative assumptions about 
ionization and metallicity corrections, the C IV absorbing 
clouds would constitute a significant baryon repository 
when compared to the galaxy stellar mass (at least ~10% 
of the M*). 

While our understanding of how starbursts affect the 
circum-galactic gas is far from complete, our new ob- 
servations show for the first time that starburst-driven 
outflows in the nearby universe can have a major impact 
on physical conditions in the CGM. This might hold the 
key to understanding whether such winds may quench 
star formation by reducing the supply of infalling cool 
gas and help us to better understand how the intergalac- 
tic medium is chemically enriched and heated. 
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Our QSO-Galaxy Sample with Information on the Foreground Galaxy and the Background QSO. 
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.473 


88, 


.7 


J120314 


.43+480316.4 


180 


.810 


48. 


,055 


0, 


,0629 


Passive 


0. 


12+0. 


12 


0.04+0.12 


10 


.78 


10 


.20 


7 


J120329. 84+480313. 5 


180. 


,87 


48. 


,05 





.817 


187. 


,4 


J122115 


.76-020009.6 


185, 


.316 


-2 


.003 


0, 


,0625 


SB/PSB 


-5, 


,06+0. 


,08 


0.89+0.09 


10 


.15 


9 


.43 


85 


J122115.32-020253.3 


185. 


,31 


-2 


.05 





.785 


197. 


,2 


J122534 


.26-025029.1 


186 


.393 


-2 


.841 


0, 


,0673 


SB/PSB 


-5, 


,73+0. 


,04 


0.88+0.04 


10 


.13 


9 


.69 





J122534.79-024757.2 


186. 


,40 


-2 


.80 





.195 


196. 


,1 


J132107 


.48+295615.3 


200 


.281 


29. 


,938 


0, 


,0723 


SB/PSB 


-3, 


,79+0. 


,11 


0.16+0.12 


10 


.48 


9 


.97 


33 


J132059. 41+295728.1 


200. 


,25 


29. 


,96 





.206 


175. 


,5 


J132150 


.89+033034.1 


200 


.462 


3, 


,509 


0. 


,0816 


Normal 


-2, 


,70+0. 


,09 


0.62+0.09 


10 


.81 


10 


.13 


16 


J132144.97+033055.7 


200. 


,44 


3, 


,52 





.269 


140. 


,2 


J140502, 


.20+470525.9 


211 


.259 


47, 


,091 


0, 


,1452 


SB/PSB 


-5, 


,02+0. 


,13 


-0.19+0.14 


10 


.43 


10 


.08 





J140505.77+470441.1 


211. 


,27 


47, 


,08 


1 


.240 


146. 


,9 


J151136 


.53+402852.6 


227 


.902 


40, 


,481 


0, 


,0945 


Passive 


0. 


04+0. 


14 


0.35+0.14 


10 


.67 


10 


.22 





J151140.33+402721.5 


227. 


,92 


40, 


,46 





.558 


177. 


,0 


J154527 


.12+484642.2 


236 


.363 


48. 


,778 


0, 


,0752 


Normal 


-2. 


,37+0. 


,19 


0.12+0.20 


10 


.50 


9, 


.87 


10 


J154530.23+484608.9 


236. 


,38 


48. 


,77 





.399 


64 


.7 


J161708, 


.92+063822.2 


244, 


.287 


6, 


,640 


0, 


,1526 


Passive 


0. 


26+0. 


11 


0.50+0.12 


11 


.39 


10 


.81 





J161711.42+063833.4 


244. 


,30 


6, 


,64 





.229 


103. 


,2 


J161913 


.50+334146.8 


244, 


.806 


33, 


,696 


0, 


,1326 


Normal 


-2, 


,23+0. 


,15 


0.24+0.16 


11 


.04 


10 


.50 


12 


J161916.54+334238.4 


244. 


,82 


33, 


,71 





.471 


150. 


,9 


J230842 


.91-091112.8 


347 


.179 


-9 


.187 


0, 


,1902 


Normal 


-2, 


,72+0. 


,11 


0.19+0.11 


11 


.09 


10 


.69 


2 


J230845.60-091123.9 


347. 


,19 


-9 


.19 





.193 


131. 


,0 



a From the latest calculation by Jarle Brinchman using the method published by Kauffmann et al. ( 2 003T) and[Salim et al. J2007T) . 

b Orientation is the angle between the QSO sightline and the major axis of the galaxy. The major axis was identify using position angle from SDSS. These values were not accompanied by 
individual errors and therefore, we estimated a mean error of 10° in orientation . 
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TABLE 2 

Table Presenting Star Formation Related Properties of Galaxies in Our Sample. 



Foreground Galaxy 


Ra Flux a 


Ha Luminosity" 


Q* 


Avg. sSFR b 


sSFRc 




(erg s _1 cm -2 ) 


(erg s _1 ) 


(s- 1 ) 


(Log (yr- 1 )) 


(Log (yr- 1 )) 


J075622 










-10.70 


J 084356 


2.3x10 x * 


7.7xlO* A 


5. ox 10°° 


-0.00 


-9.30 


J085252 


7.7xl0 -17 


1.8xl0 39 


1.3xl0 51 




-10.65 


J085254 


l.lxlO -15 


3.0xl0 39 


2.2xl0 51 




-10.85 


J085301 


3.4xl0 -16 


1.5xl0^ u 


l.lxlO 52 




-11.55 


TAArkO A A 

J092844 










-11.10 


J100801 


1.8x10 


1.2/ 10 


8.7xl0 5i 




-10.55 


J101008 


1.9xl0 -17 


3.5xl0 38 


2.6xl0 50 




-11.25 


J102846 


7.6xl0 -15 


2.5xl0 41 


1.8xl0 53 


-9.20 


-9.80 


J120314 


2.9xl0- 15 


2.8xl0 40 


2.0xl0 52 




-11.55 


J122115 


1.5xl0 -14 


1.4xl0 41 


l.OxlO 53 


-8.76 


-9.75 


J122534 


5.1xl0- 14 


5.5xl0 41 


4.1xl0 53 


-8.55 


-9.40 


J132107 


2.7xl0- 14 


3.5xl0 41 


2.5xl0 53 


-9.73 


-9.45 


J132150 


4.2xl0- 15 


7.0xl0 40 


5.1xl0 52 




-10.25 


J140502 


1.3xl0- 14 


7.2xl0 41 


5.3xl0 53 


-9.28 


-9.05 


J151136 










-11.55 


J154527 


1.8xl0- 15 


2.5xl0 40 


1.8xl0 52 




-10.50 


J161708 










-11.60 


J161913 


3.3xl0- 15 


1.5xl0 41 


l.lxlO 53 




-10.10 


J230842 


3.5xl0- 15 


3.6xl0 41 


2.6xl0 53 




-10.20 



a Extinc tion corrected flux and luminosity values as derived from MPA-JHU catalog using prescription pub- 
lished bv I Wild et all (T2007h 



Average s pecific SFR derive d based on the age of the new stellar population and burst mass fraction as 
described bv lWild etldl (120101 ). 

c Based on SDSS fiber spectra from Brinchma nn et all ([2004). These were determined from Ha and/or 
D.n,(4000) depending on availability of Ha and H/3 measurements that are uncontaminated by AGN emission. 



TABLE 3 

Lyman a and C IV Absorption Features Associated with Foreground Galaxies. 



Galaxy 


WLya res t 


a,b 


FWHM c Lya 


Log 


N L 


yot 


W C IV1548 rest 


FWHM &IV1548 


W C IV1550 rest 


Log 


Nciv 




(mA) 




(km s- 1 ) 


(Log i 


[cm 


- 2 )) 


(mA) 


(km 




(mA) 


(Log 


(cm- 2 )) 


J075622 


< 


146 


b 
















< 


255 








J084356 


437 


± 


43 


104 l +21 -° 


14.11 


± 


0.05 


< 160 






< 


215 








J085252 


279 


± 


23 


< 29.8 d 


>13.71 d 


< 199 






< 


256 








J085254 
















< 121 






< 


114 








J085301 


1452 


± 


47 


186.6+^g 3 


>15.05 e 


928 ± 86-^ 


162. 


n +56.0 
U -41.6 






14.59 


± 




J092844 


2885 


± 


132 


c-, r o + 20.'6 


>15.08 e 


779 ± 172 


129. 


+41.4 
U -31.4 


505 


± 104 


14.53 


± 


0.08 


J100801 


1694 


± 


168 


417.9+^-8 


14.69 


± 


0.02 


< 295 
















J101800 


640 


± 


84 


449.6+? 3 -j 


14.13 


± 


0.05 


< 231 






< 


221 








J102846 


1129 


± 


67 


182.8+21 4 


14.69 


± 


0.07 


1224 ± 72 


279. 


7+47.4 
' -40.5 


741 


± 92 


14.65 


± 


0.04 


J120314 


1815 


± 


129 


363.4t^;° 3 


14.79 


± 


0.03 


< 261 






< 


279 








J122115 


497 


± 


85 


143 n+ 93 - 3 


14.12 


± 


0.10 


689 ± 37 


152. 


pr+42.2 
°-33.0 


409 


± 56 


14.38 


± 


0.05 


J122534 


898 


± 


44 


1 99 9+ 21-7 


14.46 


± 


0.03 


550 ± 100 


115.* 


, + 115.8 
3 -57.9 


332 


± 61 


14.31 


± 


0.09 


J132107 


275 


± 


75 


398.6±|S:g 


13.73 


± 


0.20 


















J132150 


613 


± 


47 


161.3t£J 


14.23 


± 


0.02 


< 223 






< 


234 








J140502 
















627 ± 105 


319.* 


, + 107.6 
3 -80.5 


336 


± 57 


14.25 


± 


0.08 


J151136 


< 


127 


b 










< 248 






< 


289 








J154527 


1154 


± 


59 


281.9+ 1 1 6 5 ;° 2 


14.53 


± 


0.02 


< 134 






< 


125 








J161708 


< 


199 


b 










< 485 






< 


486 








J161913 


1847 


± 


45 


3l0.0±7;g 


14.90 


± 


0.01 








< 


94 








J230842 


66 


± 


21 


< 12.5^ 


>13.09 d 


< 154 

















a The error represented here takes into account both the statistical error as well as error in continuum fitting. 

k The limiting Lyman a equivalent width represents the 3 a error in the equivalent width over 200 km s _ 1 i.e. ±100 km s _ 1 in the rest 
frame of the foreground galaxy. For C IV we estimate limiting equivalent width around regions of ±100 km s _ 1 from the Lyman a line 
center. 

c FWHM takes into account the instrumental line spread function. 

^ Lines are unresolved and the column densities presented here were derived assuming the transitions to be optically thin. 
e The lines are saturated and the quoted column density is a lower limit. 

f The data coverage of the line is partial as it lies at the edge of the grating. Consequently, there could be much larger systematic 
uncertainty in the measurement than estimated in the error. 
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TABLE 4 

Transitions Associated with Two strongest Lyman a absorbers Associated 
with Galaxies J085301 and J092844. 



Galaxy 


Redshift 


Transition 




Wrest 


FWHM 6 




Log N 








(A) 


(mA) 


(km s -1 ) 


(Log cm 2 ) 


J085301 


0.1292 


H I Lya 


1215.67 


1452 ± 47 


187 tVr 


>15.05 c 






H I Ly 7 d 


972.54 


1023 ± 154 














Si II 


1190.42 


215 ± 23 














Si II 


1193.29 


345 ± 23 




82 ±1° 










Si II 


1260.42 


524 ± 47 






>13.8 e 






Si II 


1304.37 


165 ± 36 














Si II 


1526.71 


311 ± 63 














Si III 


1206.50 


790 ± 46 


171 +21 

1/1 -HQ 


13. 


81 ± 0.03 






Si IV 
Si IV 


1393.76 
1402.77 


423 ± 37 
240 ± 50 




92±\l 




>13.7 e 






CII 


1334.53 


633 ± 54 


1 


20 ±f 2 


14. 


68 ± 0.05 






cm d 


977.02 


1031 c ± 126 














CIV 


1548.20 


928 ± 86 


162 ±H 


14. 


59 ± 0.06 



J092844 



0.1541 



H I Lya 


1215 


67 


2987 ± 132 


H I Ly 7 d 


972 


54 


1063 ± 117 


Si II 


1190 


42 


424 ± 61 


Si II 


1193 


29 


526 ± 80 


Si II 


1260 


42 


647 ± 86 


Si II 


1304 


37 


382 ± 63 


Si II 


1526 


71 


578 ± 95 


Si III 


1206 


50 


611 ± 84 


Si IV 


1393 


76 


738 ± 77 


CII 


1334 


53 


777 ± 68 


ClII d 


977 


02 


1063 ± 118 


CIV 


1548 


20 


779 ± 172 


CIV 


1550 


77 


505 ± 104 



515 



+ 21 
-20 



73 



+ 8 
-7 



82 
243 
106 



+42 
-28 
+46 
-39 
+ 24 
-20 



126 



+41 
-31 



>15.08 c 



>14.13 e 



>13.45 e 
14.03 ± 0.05 
>14.56 e 

14.53 ± 0.08 



a The error takes into account both the statistical error as well as error in continuum fitting. 
b The FWHMs were obtained after correcting for the appropriate line spread function of grating 
G140L. 

c Saturated lines. 

^ The transition lies in the Segment B of the G140L grating and wavelength calibration is 
highly uncertain in this part of the spectrum. Therefore, we do not use these transition for 
fitting Voigt profile. 

e Lines are unresolved and the column densities presented here were derived assuming the 
transitions to be optically thin. 
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TABLE 5 

Absorption-line Measurements from Stacks for the Three Sub-samples. 



Galaxy 


Transition 


(A) 


w r a i t 

(mA) 


FWHM c 
(km s- 1 ) 


Log N d 
(Log (cm- 2 )) 


SB/PSB 


H I Lya 
CII 
CIV 
Si III 
Si IV 


1215.67 
1334.53 
1548.20 
1206.50 
1393.76 


590.4 ± 30.0 

<60.1 
559.6 ± 42.3 
125.2 ± 19.5 

<44.0 


171 t\l 

278 +^ 
205 ±11 


>14.19 e db 0.03 
<13.5 
14.20 ± 0.04 
12.85 ± 0.12 
<12.71 



Control H I Lya 1215.67 881.3 ± 21.5 331 ±{l >14.32 e ± 0.01 

(Normal) CII 1334.53 <40.1 <13.32 

CIV 1548.20 <39.5 <13.00 

Si III 1206.50 255.1 db 18.7 332 13.10 db 0.05 

Si IV 1393.76 111.2 ± 20.0 146 ±l% 8 13 - 12 ± °- 10 



Control 


H I Lya 


1215.67 


946.2 db 31.7 


425 


+ 23 
-22 


>14.33 e ± 0.02 


(Passive) 


CII 


1334.53 


213.5 db 31.7 


383 


+88 
-65 


14.06 ± 0.08 




CIV 


1548.20 


487.5 db 56.0 


282 


+ 70 
-57 


14.13 ± 0.06 




Si III 


1206.50 


264.0 db 28.7 


274 


+83 
-63 


13.12 db 0.06 




Si IV 


1393.76 


<167.7 






<13.30 



a The lower limits were derived using data within ztlOOkm s 1 of the expected line 
center. 



The error represented here takes into account both the statistical error as well as 
error in continuum fitting. 

c The FWHMs were obtained after correcting for the appropriate line spread function 
of grating G140L. 

d Column densities for the non detections were derived by assuming a Doppler 
parameter, b = 60km s — 1 . 

e Some of the absorption features that were included for making the stacks were satu- 
rated. Therefore, the column density represented here is a lower limit. 



